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Summary 

A technique is presented whereby the mutual inter- 
action between a small number of elements in a planar 
array can be interpolated and extrapolated to accurately 
predict the combined interactions in a much larger array 
of many elements. An approximate series expression is 
developed, based upon knowledge of the analytical char- 
acteristic behavior of the mutual admittance between 
small aperture antenna elements in a conducting ground 
plane. This expression is utilized to analytically extend 
known values for a few spacings and orientations to other 
element configurations, thus eliminating the need to 
numerically integrate a large number of highly oscillat- 
ing and slowly converging functions. This paper shows 
that the technique can predict very accurately the mutual 
coupling between elements in a very large planar array 
with a knowledge of the self-admittance of an isolated 
element and the coupling between only two elements 
arranged in eight different pair combinations. These 
eight pair combinations do not necessarily have to corre- 
spond to pairs in the large array, although all of the indi- 
vidual elements must be identical. 

Introduction 

The electromagnetic performance (e.g., impedance, 
radiation pattern, and scanning characteristics) of an 
array antenna sometimes is influenced significantly by 
the combined mutual interaction of the local fields of the 
array elements. A good antenna design would include an 
accurate account of these interaction effects. 

The radiation characteristics of an array antenna can 
be obtained from a knowledge of the properties of a 
single element, the array geometry, and the combined 
interaction of all elements of the array. The combined 
mutual interaction of all elements can be obtained from 
the scattering matrix in which the size of the matrix is 
determined by the number of elements in the antenna 
array (e.g., an antenna containing N elements would be 
represented by an N x N scattering matrix). The complex 
coefficients of the scattering matrix can be obtained from 
the complex admittances of all element pair combina- 
tions of the antenna array, either by calculations for ele- 
ments that can be adequately modeled or by 
measurements for more complicated elements. 

Since these interactions are typically a near-field 
phenomenon, a theoretical formulation will usually 
involve the numerical evaluation of quite complicated 
integrations, resulting in significant computational time 
for an array of many elements. The measurement of 
mutual interactions between all pair combinations within 
a large array can also become very laborious and time 
consuming. A technique is presented that utilizes the 


characteristic analytical behavior of the coupling 
between two elements, in conjunction with known values 
of coupling for a few configurations, to accurately pre- 
dict the mutual coupling for any configuration of these 
elements in an array. 

Method 

The mutual admittance between two copolarized 
apertures in a common, flat, perfectly conducting ground 
plane is postulated as 



( 1 ) 


The geometry of two rectangular apertures is shown in 
figure 1, where R is the center-to-center spacing of the 
apertures, <(> is the relative angular orientation of the aper- 
tures, and k = 2 n/X, with X being the operating wave- 
length in free space. The polarization direction of the 
apertures is assumed such that <|> = 0 represents coupling 
in the magnetic field direction (i.e., //-plane) and 0 = k/2 
represents coupling in the electric field direction (i.e., 
E-plane). Aj,A 2 , ■•■Ms are complex coefficients to be 
determined from known values of F 12 for pairs of aper- 
tures with known relative positions. 

The reasoning that influenced the postulated form of 
equation (1) follows: 

(a) Previous observations have shown that the mutual 
coupling between widely spaced apertures in a 
common ground plane varies inversely as the sepa- 
ration distance in the E-plane and inversely as the 
square of the separation distance in the //-plane, 
with a phase variation that is linear with separation. 

(b) The fields near a radiating source vary inversely as 
the cube of the distance away from the source. 

(c) For a fixed separation, the mutual coupling varies 
in a sinusoidal manner with relative angular 
position. 

These basic analytical properties are also illustrated 
by the closed-form expression for mutual admittance 
between dominant-mode circular apertures (ref. 1). It 
is anticipated that the mutual admittance between any 
other pair of apertures would exhibit similar character- 
istics versus separation and angular position; thus the 



postulated expression in equation (1) is expected to apply 
to a wide variety of elements in planar arrays. However, 
the coefficients in equation (1) must be determined for 
each unique type of element pair. 

In postulating the form of equation (1), the sinusoi- 
dal functions were chosen such that the //-plane and 
Zs-plane mutual-coupling values could be determined 
independently. If one knew the complex mutual admit- 
tance between two elements for two spacings in the 
//-plane (<() = 0), the coefficients A j and A 2 could be 
determined. Likewise a knowledge of the mutual admit- 
tance for three spacings in the E-plane (<j) = tc/ 2) would 
be adequate for determining the coefficients A 3 , A 4 , 
andA 5 . The first five terms are usually sufficient for 
approximating the mutual admittance for other element 
pair combinations because of the cylindrical interpola- 
tion for a specified value of R by the sin 2 <|) and cos 2 <|) fac- 
tors. The inclusion of three additional terms (with 
multiplying factor sin 2 2<)>) yields even more accurate 
results for values of 0 < 0 < n!2. Eight distinct element- 
pair combinations are used to establish a set of eight 
independent equations to be solved for the eight 
unknown coefficients by matrix inversion. These pair 
combinations do not necessarily have to include the 
//-plane and £-plane, although these would be logical 
choices. 

In obtaining sample data for the evaluation of coeffi- 
cients in equation (1), some rather loose (although 
logical) restraints on the selection of element-pair combi- 
nations should be observed in order to obtain realistic 
and accurate interpolation/extrapolation values. No more 
than three element spacings for the same value of <(> 
should be used for sample data (and no more than two 
sample spacings for <J) = 0). Sample data for very close 
spacing (where the dominant terms are 1 /R 3 ), sample 
data for large spacing (where the dominant terms 
are UR ), and sample data for intermediate spacings 
should be used. 

The coefficients in equation (1) could be determined 
by either calculations or measurements. Calculations are 
appropriate for elements that can be adequately modeled 
analytically, such as waveguide-fed apertures or cavity- 
backed slots. However, more complicated elements do 
not lend themselves easily to straightforward analyses; 
therefore, a measurement approach may be more appro- 
priate in some cases. 

If measurements are used to determine the coeffi- 
cients, the simplest approach is to measure the reflection 
coefficient and coupling coefficient (i.e., scattering 
matrix) for a pair of elements (i.e., a two-element array) 
for various spacings and positions. The mutual admit- 


tances can then be determined from the relationship 
between the admittance and scattering matrices. 

\s] = fm-iymm + mr 1 (2a) 

[y] = [[/]-[$]][[/] + [S]]- 1 (2b) 

where [5] is the complex scattering matrix, [/] is the 
identity matrix, [y] is the complex normalized admit- 
tance matrix, and [ J _1 indicates matrix inversion. The 
coefficients of the admittance matrix in equations (2a) 
and (2b) are normalized to the ith individual modal- 
characteristic admittance for a waveguide of the same 
cross-sectional dimensions as the ith aperture, = YgfYj, 
where Y i is the modal-characteristic admittance for the 
ith aperture and where subscripts i and j refer to the ith 
row and yth column of the matrix. The measured mutual 
admittance between two identical apertures can then be 
determined from 

Y \2 - 25 12 

y t ~ y 12 " (i +s n ) 2 -s2 2 (3) 

where S n and S 12 are the measured complex reflection 
and coupling coefficients, respectively, of the two- 
element array. 

Results 

For purposes of this paper, calculations based upon 
the numerical evaluation of integral expressions for cir- 
cular, rectangular, and square apertures are used to deter- 
mine the coefficients in equation (1) and to evaluate the 
accuracy of the utility of equation (1). The integral 
expressions for mutual admittance of circular and rectan- 
gular apertures can be obtained from the general formu- 
lation in reference 2. The computer program used for 
circular aperture results is documented in reference 3, 
and the computer program used for rectangular and 
square apertures is documented in reference 4. 

Figures 2 and 3 show the arrangements of elements 
used for determination of the complex coefficients in 
equation (1) for rectangular and circular apertures. The 
geometries for figures 2 and 3 were selected to satisfy 
the loose constraints mentioned earlier. These arrange- 
ments were found to give very accurate interpolated/ 
extrapolated results. It was also found that these specific 
arrangements were not critical to obtaining accurate 
results. Some reasonable deviations from these geome- 
tries yielded insignificant deviations in the final results, 
provided that the loose constraints were met. The mutual 
admittance between the element at the origin and each of 
the other eight elements was calculated from the respec- 
tive computer programs, RWG and CWG, by numerical 
integration. These calculated values were used to solve 
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the corresponding set of eight simultaneous equations for 
the coefficients. These coefficients were then used in 
equation (1) to predict the mutual admittance between 
two elements for a wide range of spacings and orienta- 
tions, and comparisons were made with corresponding 
results obtained by direct numerical integration. 

Figures 4 through 9 show a comparison between 
rectangular aperture mutual admittance obtained by 
equation (1) and the values obtained by numerical 
integration. The sample data points are those used for 
interpolation/extrapolation. These data show both the 
real and imaginary parts of the complex mutual admit- 
tance between two identical apertures as a function of the 
aperture spacing at specific orientation angles <f). Fig- 
ures 10 through 14 show the results for circular aper- 
tures. It is readily evident that very accurate results can 
be obtained for any orientation and spacing of apertures 
with eight sample data points. 

Interpolated/extrapolated results obtained from five 
sample data points (//-plane and £-plane only and coeffi- 
cients A | through A 5 in eq. (1)) for circular apertures 
were indistinguishable from those in figures 10 through 
14. A very slight, although still insignificant, deviation 
was observed in the <f> = tc/ 4 plane results when using 
only the five sample data points for rectangular aper- 
tures; therefore, use of all eight coefficients is recom- 
mended if sufficient sample data are available. 

Figures 1 5 through 1 8 show a comparison between 
interpolated/extrapolated mutual-coupling results with 
those obtained by direct numerical integration. Very 
accurate results are indicated for extrapolation to spac- 
ings of 15 wavelengths and greater and for very low cou- 
pling levels, although the sample data spacing was 
limited to about 5 wavelengths or less, implying that the 
technique could be used for very large arrays of many 
elements. The results using direct numerical integration 
for the square apertures exhibits an oscillatory behavior 
for wider spacings, which indicates possible errors in the 
calculations. This accurate extrapolation to large spac- 
ings was indeed verified in the expanded plots in fig- 
ures 19 through 22, where the computer codes were mod- 
ified for improved accuracy in the numerical integration, 
with a corresponding increase in computation time. With 
improved numerical accuracy, the amplitude of the oscil- 
lations became smaller and closer to the extrapolated 
data. These results illustrate an important aspect of the 
present technique. As the spacing between array ele- 
ments increases, the numerical-integration approach 


becomes less accurate and more computationally inten- 
sive; however, the algebraic extrapolation technique 
maintains extreme accuracy without increasing computa- 
tion time. 

Concluding Remarks 

An algebraic interpolation/extrapolation technique 
for predicting the mutual interaction between elements of 
a large planar phased-array antenna has been presented. 
The technique has been shown to be extremely accurate 
although it only uses eight sample data points. The 
technique yielded very accurate results for large ele- 
ment spacings and low coupling levels, as opposed to 
numerical-integration techniques that exhibit oscillatory 
errors and large increases in computation time. Although 
the technique has been verified only for dominant- 
mode circular and rectangular aperture elements and 
the interpolation/extrapolation equation may take on a 
slightly different form for other elements or apertures 
with different modes present, this general approach 
should be applicable to a wide variety of array element 
types. The present extrapolation technique may also pro- 
vide significant advantages over methods of direct mea- 
surement of coupling between widely spaced elements 
because extreme care must be exercised in order to 
obtain accurate measured results for low coupling levels. 
An attractive feature of the present technique is that it 
utilizes and extrapolates from data for strongly coupled 
elements that are closely spaced, the very situation where 
measurement methods or numerical integration methods 
can provide it with accurate sample data. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
May 21, 1996 
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Figure 2. Rectangular element arrangement for determination of admittance coefficients. 
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Figure 5. Mutual admittance between two rectangular apertures (0.8A. x 0.4A). <)) = n/8. 



Figure 6. Mutual admittance between two rectangular apertures (0.8 A. x 0.4A). arctant)) = 1/2. 
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Figure 8. Mutual admittance between two rectangular apertures (0.8A, x 0.4A,). <() = 3 ji/8. 
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Figure 14. Mutual admittance between two circular apertures (0.6A. diameter). <|) = rc/2; £-plane. 
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Figure 15. Mutual coupling between two circular apertures (0.6^ diameter). 
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Figure 16. Mutual coupling between two circular apertures (0.8^ diameter). 
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Figure 18. Mutual coupling between two square apertures (0.6A, x 0.6A,). 
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Figure 19. Amplitude of //-plane coupling between two square apertures. 



Center- to -center spacing, X 

Figure 20. Phase of //-plane coupling between two square apertures. 
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Figure 21 . Amplitude of £-plane coupling between two square apertures. 
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Figure 22. Phase of £-plane coupling between two square apertures. 
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